ON ONE-POINT FUNCTIONS OF DESCENDANTS IN 
SINE-GORDON MODEL 



M. JIMBO, T. MIWA AND F. SMIRNOV 



Abstract. We apply the fermionic description of CFT developed in our previous 
work to the computation of the one-point functions of the descendent fields in the 
sine-Gordon model. 



1. Introduction 

The sine-Gordon (sG) model is the most famous example of two-dimensional inte- 
grable Quantum Field Theory (QFT). The sG model is defined in two-dimensional 
Minkowski space with coordinate x = {xo,Xi) by the action 

(1) ^sG= / (-^(a,C^(x))2 + -4^COs(Mx))|rf2x. 

J I Ion sm^Trp"') J 

The normalisation of the dimensional coupling constant in front of cos(/3£^(x)) is 
chosen for future convenience. This model has been a subject of intensive study 
during the last 30 years. First, by semi-classical methods the exact spectrum was 
computed, the factorisation of scattering was predicted, and the exact S'-matrix was 
found for certain values of the coupling constant (in the absence of reflection of 
solitons) [2113]. 

The most significant further results were found by the bootstrap method. In [U [S] 
the exact S-matrix was found, and in [HI E] the exact form factors were computed 
for the energy-momentum tensor, topological current and the operators e^^^'^^^\ 
e='=^v'(x)_ Then the latter result was generalised to the operator e^"''^^^'^ with an 
arbitrary a [8j, using the methods which go back to the study of the closely related 
XXZ spin chain [9, lOj. 

It should be said that many important technical and conceptual methods of the 
modern theory of quantum integrable models originate in the quantum inverse scat- 
tering method [HI [T2] . It provides a clear mathematical interpretation of the work 
by R. Baxter [13] . In particular, in [12\ the scattering matrix of the sG solitons was 
reproduced using the Bethe Ansatz. 

The knowledge of form factors allows us to write a series representation for the 
two-point function 

(2) ^giai^(x)gm2VP(0)^^^ _ 

In this paper we shall consider only the space-like region x^ < 0. We shall use a lat- 
tice regularisation which breaks the Lorentz invariance. So we shall take x = {0,x), 
and use the notation ip{x) = ip{0,x). In this case the integrals over the form factors 
are rapidly convergent. It is rather hard to give a mathematically rigorous proof 
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of the convergence of the series, but nobody doubts that they converge. Actually 
the convergence was proved [2] for a certain particular reduction of the sG model, 
known as the scaling Lee- Yang theory. So, the situation looks completely satisfac- 
tory. However, the series over the form factors converge slowly in the ultraviolet 
region (for small values of — x^). To give an efficient description to its ultraviolet 
behaviour remained a problem largely open for some time. Before describing the 
way of solving this problem let us explain how it is solved in a particular case. 

It is known that at = 1/2 the sine-Gordon model is equivalent to the free the- 
ory of a Dirac fermion. The correlation functions of ([2]) are non-trivial. The sG field 
^p{x) is bilinear in fermions, and one has to compute the correlation function of two 
exponentials of bilinear forms. The result is obtained in [TB] , generalising the sem- 
inal work [TB] on the scaling Ising model. Namely, it is shown that the correlation 
function satisfies an equivalent of a Painleve III equation. The form factors are very 
simple in this case, and the form factor series coincides with the Fredholm determi- 
nant representation of this solution. Still the problem of describing the ultra-violet 
behaviour is non-trivial. It amounts to finding the connection coefficients for the 
Painleve equation, which is done by studying the Riemann-Hilbert problem pTTj. 
Through this analysis, one draws an important conclusion. Since CFT completely 
describes the ultra-violet limit, one might naively expect that the asymptotic be- 
haviour of the two-point function may be obtained via the perturbation theory. Such 
an assumption would imply that the dependence on the mass scale is analytic in /x^ 
which has the dimension (mass)^*^^~^ However this is not the case even for the 
Painleve solution. 

Correct understanding of the conformal perturbation theory is one of the most 
important problems in the theory of quantum integrable systems. This problem was 
studied in [18j. In fact, the naive perturbation theory suffers from both ultraviolet 
and infrared divergencies. The idea of ^8J is to absorb all these divergencies into 
non-perturbative data: one-point functions of primary fields and their descendants. 
Once it is done, the remaining task is a convergent version of conformal perturbation 
theory. So the problem is divided into two steps. The first one requires some non- 
perturbative information. The second one is of genuinely CFT origin: actually, it 
is reduced to the computation of some Dotsenko-Fateev Coulomb gas integrals with 
screenings P^ . 

In principle, the procedure described in [TS] provides an asymptotic series in the 
ultra-violet domain which agrees with the structure expected from the Painleve 
case: it is not just a power series in /x^, but includes non- analytic contributions with 
fractional powers of . 

So the main problem is to compute the one-point functions. The first important 
result in this direction was achieved in [20], where an exact formula for the one- 
point functions of the primary fields was conjectured. Then by several ingenious 
tricks (such as going form sine-Gordon to sinh-Gordon and back) a procedure was 
described in [211 [221 [23] , which must in principle allow us to compute the correlation 
functions of the descendants. Unfortunately, this procedure involves certain matrix 
Riemann-Hilbert problem which has not been solved in general up to now. At the 
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same time, this way of computation looks very indirect, and involves steps which are 
hard to justify. Still, the predicted result for the first not-trivial descendant [23] is 
quite remarkable. Even though it was obtained by a complicated and non-rigorous 
procedure, it was checked against many particular cases. So we do not doubt in its 
validity. It will be used to check the result of the present paper. 

Let us mention here a deep relation between the sG model and the $i,3-perturbation 
of c < 1 models of CFT [21] . On the formal level the relation is simple. In the action 
([T]) one can split cos(/3(^(x)) = |(e~*'''^*^'') -|- e'^^'^^^^) and consider the first term as a 
part of the Liouville action and the second as a perturbation. The Liouville model 
with an imaginary exponent is nothing but CFT with the central charge 

c=l-6(/3-l//3)2, 

and e^^'^^^^ is the field <l>i 3(x). So, formally there is no difference between the sG 
model and the $i^3-perturbation of c < 1 CFT. The situation becomes interesting 
in the case of rational It was shown in [251 that in the computation of the 
correlation functions of e"^^^^^ {n = 1, 2, 3, ■ ■ ■ ), a restriction of degrees of freedom 
takes place for the intermediate states of solitons. The mechanism of reduction is 
similar to the RSOS restriction for vertex models. This phenomenon is known as 
the restriction of the sG model. For example, if = 3/4 the solitons reduce to 
Majorana fermions, and the restricted model is nothing but the scaling Ising model. 

In the series of papers, [271 ISH [291 [I] which we will refer to as I, II, III, IV, re- 
spectively, we studied the hidden fermionic structure of the XXZ spin chain. In 
particular, in IV the relation between our fermions in the scaling limit and c < 1 
CFT was established. The long distance behaviour of the XXZ model and the short 
distance behaviour of the sG model are described by the same CFT: free bosons with 
the compactification radius For the XXZ model, we use the coupling parameter 
V related to via 

This identification is used when relating the results of IV to those of [30l [31] , which 
were important for us. The relevant CFT has the central charge 

c= 1-6 . 

1 - V 

Here we do not consider the peculiarity of a rational z/, so using the c < 1 CFT means 
just the usual modification of the energy-momentum tensor. In the ultraviolet limit, 
the sG model is described by two chiral copies of CFT. We use the notation 

= e T^4f v'W} . 

The field il3ifi{x) splits into two chiral fields 2(/?(a;) -|- 2ip{x). (We follow the normal- 
isation of the fields ip{x), (^(x) given in ^Ol [31] and in IV.) Our goal is to compute 
the one point functions 

(P({i_™})P({L„})$„(o)),G 



(<fa(0)).G 
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where P( {!_,„}), -P({l_m}) are polynomials in the generators of two chiral copies of 
the Virasoro algebra with the central charge c. 

The universal enveloping algebra of the Virasoro algebra contains the local inte- 
grals of motion i2k-i (ii = l-i) [32] which survive the $i^3-perturbation. Clearly, 
the one-point functions of descendants created by them vanish. We assume that it 
is possible to write any element of the chiral Verma module generated by $o(0) as 
Pi({i2A:-i})-P2({l-2fc})*&a(0). So, actually we are interested in computing 

(P({l-2m})P({l-2m})$»(0)).G 
($a(0)).G 

For the simplest non-trivial case l_2l_2$a(0) the problem was solved in [23]. If one 
wants to consider the descendants created by the Heisenberg algebra it is easy to do 
using the formulae 

(1 - u)T{z) =: ip'izy : +Mz), (1 - u)f{z) =: ^'(z)^ : +u^{z) . 

In this paper we shall consider only the domain < a < 2, but the final results 
allow analytical continuation for all values of a. 

The whole idea of IV is that the usual basis of the Verma module is not suitable 
for the perturbation, and we have to introduce another one. For the moment we 
consider only one chirality. Working modulo action of the local integrals of motion 
the new basis is provided by uncharged products of two fermions 02k-i a-iid 'y^^^i ■ 
The fermions respect the Virasoro grading: 

ri /qCFT*_,CFT*i (o^^Oh /qCFT*^,CFT* 

[lo, P2j-i l2k-i J = (2j + 2/e - 2)/^2i-i 72fc-i • 

We have 

(3) /3?r 7?-^^*$.(0) = {P,T?- ({1-2.}) + $«(0) , 

where /"*" and /~ are ordered multi-indices: for / = (2ri — l,---2r„ — 1) with 
ri < ■ ■ ■ < r„, we set 

/oCFT* _ AjCFT* /oCFT* ^,CFT* _ _,CFT* _,CFT* 
Pi — P2ri-1 ■ ■ ■ P2r„-1 > ~ ' 2r,j-l ' ' ' T2ri-1 i 

and similarly for /3j-, 7^-. We require = In the right hand side of ([3]), 

P|™^_({1_2A:}) and P^+^_ ({l_2fc}) are homogeneous polynomials, the constant da is 
given by 

da = "^^^^^{a - 1) = |V(25-c)(24A„ + l-c), 
and the separation into even and odd parts is determined by the reflection (3'^+'^*^^^'^* O 

aCFT* CFT* 

Pi- I1+ ■ 

The coefficients of the polynomials Pfl^^. and Pj+*^- are rational functions of c 
and 

a{a — 2)u'^ 
4(1-.) 

only. The denominators f actor ise into multipliers Aq, + 2k, k = 0,1,2,---. Exact 
formulae up to the level 6 can be found in IV, Eq. (12.4). In particular, on the level 
2 we have P^^^'^^ = 1, P°f^ = 0. The transformation ([3]) is invertible. 
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The operators /S^k-i and 7^^^^^ were found as the scahng hmit of the fermions 
which create the quasi-local fields for the XXZ spin chain. This became possible after 
the computation of the expectation values of the quasi-local fields on the cylinder 
(see III). Actually, two different scaling limits are possible, and the second one 
provides the fermionic operators for the second chirality: ^2!-* ^2k-i- The 
same story repeats for these operators, in particular, we have 

(4) 7F-^^*$«(0) = {P|rf- ({1-2.}) - d^P^fj4{Uk})} $.(0) . 

The main statement of this paper is that in the fermionic basis the sG one-point 
functions are simple: 

//□CFT* r;FT*/qCFT* CFT*^ fn\\ 



(5) 



($a(0))sG 

^ V 2i/ V 2n-lG/+ 2n-le/- 



Here 



Gr,{a) = {-ir-\in-l)\) 



+ ^(2n-l))r(l-f-^(2n-l)) 
r(l-f-^(2n-l))r(f + ^(2n-l)) 



|/| stands for the sum of elements of /, and M is the mass of soliton. We recall 
that #(/^) = #(/^) is required in order to stay in the same Verma module, and 
#(/+) = #(/-) follows. 

Using ([5]) we can find all the one-point functions of descendants. Up to the level 6 
the results of IV can be used. In particular, at the level 2 we find a perfect agreement 
with the formula (1.8) of [23] after the identification: i] = a — 1, ^ = To proceed 
to levels higher than 6 one should perform further computations in the spirit of IV. 

Let us explain how we proceed in justification of the main formula For CFT 
we use the lattice regularisation by the six vertex model (equivalently XXZ spin 
chain). For this model we use the fermionic description of the space of quasi- local 
operators found in I, II. On the lattice we have creation operators b*(C),c*(C) and 
annihilation operators h{(), c{(). The most honest way to proceed to the sG model 
would be to regularise it via the eight vertex model, and then to consider the scaling 
limit. On this way we would meet two problems. The first is the U{1) symmetry 
which is broken in the eight vertex model. It is hard to introduce the lattice analogue 
of $Q, with arbitrary a. This difficulty may be overcome by going to the SOS model. 
The second problem is conceptually more difficult. For the elliptic i?-matrices we do 
not have an analogue of the g-oscillators which is crucial for the construction of our 
fermionic operators. For the moment we do not know how to attack this problem. 
Let us notice, however, that in the case a = an analogue of bilinear combinations 
of the annihilation operators b((^)c(^) exists. It is defined in the papers [33t [M]. 

So, having problems with the eight vertex model we are forced to take another 
approach. In the paper [35] the sG or the massive Thirring models was obtained 
as a limit of the inhomogeneous six vertex model with the inhomogeneity Co (see 
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section [3] below). Notice that, for this hmit to make sense, one has first to consider 
the finite lattice on n sites, and then take the limit (^o ^ oo, n — )■ oo in a concerted 
way in order that the finite mass scale appears. But exactly this kind of procedure 
became very natural for us after we had computed in III the expectation values 
of quasi-local operators on the cylinder. The compact direction on the cylinder is 
called the Matsubara direction, and its size n is what is needed for considering the 
limit in the spirit of [35] . 

The paper is organised as follows. In Section |2] we review our previous paper IV, 
and explain how to obtain the fermionic description for two chiral CFT models from 
the XXZ spin chain (six vertex model). In SectionOwe introduce the inhomogeneous 
six vertex model and consider the continuous limit which produces the sG model 
according to [33] . We derive the one-point functions using the fermionic description 
of ultra-violet CFT. 



2. Two SCALING LIMITS OF THE XXZ MODEL AND TWO CHIRALITIES 

Our study of the XXZ model is based on the fermionic operators defined in I, II. 
This definition allowed us to compute in III the following expectation value. Con- 
sider a homogeneous six vertex model on an infinite cylinder. Let Tg m be the mon- 
odromy matrix, where S refers to the infinite direction (called the space direction), 
and M refers to the compact circular direction (called the Matsubara direction). 
We use n to denote the length of the latter. We follow the notations in IV. For a 
quasi- local operator g2o5(o)Q ^^le spacial lattice, we consider 



(6) Z^{g2"5(o)o| 



TrsTrM(Ts,Mg2-^+2.5(0) 

The generalisation of this functional Z^'* was introduced in IV. In the scaling limit, 
the introduction of s in this functional amounts to changing (screening) the back- 
ground charge at a; = — oo by —2s^^. It enables us to deal with the special case 
of the functional for which the effective action of local integrals of motion becomes 
trivial. 

The quasi-local operators are created from the primary field q^""^^^^ by action of 
the creation operators t*(^), b*(^), c*{Q. Actually they act on the space 

oo 
s=—oo 

where "\Na-s,s denotes the space of quasi-local operators of spin s with tail a — s. 

In this paper we shall consider the subspace W[°|j^ of the space W(") which are 
created from the primary fields only by fermions b*, c* (see ( IT2|) below). On w[°^j^ 
the functional Z^'* allows the determinant form which is convenient to summarise 
as 



{(1 ^1 Trs(g2-^(0)) 
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where 

Iff dC^ d^"^ 

^" " J2^ J J ^-t,n(C,0c(0b(C)^ 
r r 

where the contour F goes around = 1. 

The function Wrat.nlC^O is defined by the Matsubara data (see III, [37], IV). 
Besides the length of the Matsubara chain n it depends on the parameters k, a, 
s and on possible inhomogenieties in the Matsubara chain. However, we mark 
explicitly only the dependence on n which is the most important for us here. 

We have ignored the descendants created by t*(C). Actually, they do not give any 
non-trivial contributions in the limit n — )■ oo, which we shall be interested in. Then 
is automatically reduced to the the quotient space: 

= w("V(t*(C) - 2)W(") , 

which is obviously isomorphic to y^f^rm ^ linear space. 

We shall not repeat the definition of cJrat.nlC^O because in the present paper 
we shall use only very limited information about it. Let us explain, however, the 
suffix "rat": the function (iC/C)"'^rat,n(C) is a rational function of and The 
really interesting situation occurs when n — )■ cxd. In that case the Bethe roots for the 
transfer matrices in the Matsubara direction become dense on the half axis > 0. If 
we do not introduce additional rescahng as described below, the function uJra.t,n{CyO 
goes to the simple limit: 

(7) a;rat,„(C, 4a;o(C/e, a) + Va;(C/e, «) , 

where 

ioo 

, . f sin ^((1 — z/)m — a) 

8 4u;o C,« =- / C-^r 

J sm |(u — a) cos 

—ioo 

Vu{C, a) = -^{(q, a) + ij{(q-\ a) + 2iC tan(^) , 



2{e - 1) 



The reason for extracting the elementary Vw is due to the fact that the function ujq 
satisfies the relation typical for CFT 

(9) u^{C.a) = Uo{C\2-a). 

Notice that (^/C)°i^o(C/^) ct) is not a single-valued function of and So, the 
property of rationality is lost in the limit. We define VLq and in the same way 
as fin replacing Wrat,n(C;0 by respectively AujQ{C,/i,a) and Va;(C/^,«)- 

Following IV we denote the original creation operators introduced in II by b*^^(C) 
and c*^t{C)- They satisfy the property: 

Trs (b:,,(C)(X)) = 0, Trs «,,(C)(X)) = , 
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for all quasi-local operators X. In the present paper it is useful to replace these 
operators by the following Bogolubov transformed ones: 

(10) b*(C) = e-^%:,,(C)e^^ c*(C) = e-^^<,(C)e^^ . 

Obviously, the functional Z^'* calculated on the descendants generated by these op- 
erators is expressed as determinant constructed from function (jJrat,n(C) O^^'^lC/^i ct) 
which in the limit n — >■ cxd goes to 4a;o(C/^, ct)- 

As it is explained in IV, the operators C~°'^o{0 ^"^^ ("'^oiO rational functions 
of as far as they are considered in the functional ^n'*, and they have the following 
behaviour at = 0: 

oo oo 
bo(C) = C bsj.j,(jpjjj, Cq((^) = C ^ ''^screenj • 

For b*j,j,(,gjjj, c*j.j.gg^ .^we have used the suffix "screen" to stand for "screening" in view 
of its similarity to the lattice screening operators used in IV. 

In IV, another set of operators b*(C), c*{() is obtained from h*^^{(), c*^^{() by a 
kind of Bogolubov transformation which contains t*((^). We shall not write explicitly 
this Bogolubov transformation, but only the one relating b*(C), c*{() to bQ(^), Cq(^), 
both acting on the quotient space Wquo because only these operators are used in 
this paper. The operators b*(C), c*{() are important because Z^* vanishes on their 
descendants. Notice that we do not allow k, s to grow together with n. In that case 
the dependence on k,s disappears for n = oo. 

Remark. The defining equations for uJrat,n{C,0 given in IV imply that the limit 
(171) has a very general nature. Namely, the result is independent not only of k, s 
but also of inhomogeneities in the Matsubara chain. The situation is similar to that 
for the ^-matrix which is the same for homogeneous, inhomogeneous XXZ chains 
or even for the sG model. 

If we consider the Bogolubov transformation which connects the operators b*((^), 
c*{() and bo(C), Cq(C) as acting on the quotient space Wquo, it reduces to 

(11) b*(0 = e-^°b*(C)e^«, c*(C) = e-^°c*(C)e^«. 

We catch operators acting on W*^"^ by developing b*((^), c*{() and b(C), c{Q 
around the point ('^ = 1: 

oo oo 

(12) b*(c) ^ Y.^e - irx, c*(c) ^ Y.^e - ly-'c; , 

p=i p=i 

oo oo 

(13) b(c) = - <o = - 1)"% • 

p=0 p=0 

These operators are (quasi-)local in the sense of II, Section 3.3, while b*(,j,ggjjj and 
<creenj ^rc highly "non-local" . 

The operators b*, c* create quasi-local operators in the sense of II, Section 3.1 by 
acting on the primary field. Here we slightly change the notation compared to II, III, 
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where the Fourier coefficients are defined after removing an overall power (^^". The 
reason was that we wanted that the result of action of b*, c* is a rational function 
of g, g". This rationality property is irrelevant in this paper, and extracting 
may even cause a confusion. 

So far we have been discussing the simple limit n ^ oo. Now we discuss the 
scaling limit to CFT. In IV we studied the scaling limit of the homogeneous XXZ 
chain on the cylinder. The key idea is to consider first of all the scaling limit in the 
Matsubara direction. Namely, denoting the length of the Matsubara chain by n and 
introducing the step of the lattice a we consider the limit 

(14) n — 7- oo, a — 7- 0, na = 2itR fixed . 

The requirement is that if we rescale the spectral parameter as ( = (Ca)^A, the 
Bethe roots for the transfer matrix in the Matsubara direction which are close to 
= remain finite. The constant C is chosen to have an agreement with CFT as 

(15) c= /Jj^. m^- 

Next we consider the scaling limit in the space direction. We conjecture that 
under the presence of the background charges effected by the screening operators, 
the lattice operator g2KS-f2aS'(o) g^gg ^.^ ^^le limit $i_K'(~^)'^a(0)$i+K(oo) where 

/ 1 - 
K = K + a + 2s . 

V 

Furthermore, we have shown that in the weak sense the following limits exist 

(16) 2/3*(A) = lim h\{CaY\), 27*(A) = lim c\{CaY\) . 

a— J-O a— s-0 

The operators /3*(A), 7* (A) have the asymptotics at A^ — i- 00: 

00 00 

(17) /3*(A) ^ ^'"^ (31^-1^ 7*(A) ^ ^'"^^h-i ' 

i=i i=i 

where /32j_i) 72j-i ^ct between different Verma modules as follows 

1*23-1 ■ ^a+2l^{s+l)®^c. V„+2i^, ®V«. 

The action on the second component is trivial. The functional Z^'^ turns in this limit 
into the three-point function for c < 1 CFT as explained in IV. The identification of 
descendants created by (3* and 7* with Virasoro descendants is made by studying 
the function 

ujR{X,fi) = h lim Wrat,n((Ca)''A, (Ca)» - Va;(A//i, a)) , 

4 scaling 

where limgcaUng refers to the scaling limit ( IT4|) . Contrary to the simple-minded limit 
n — )■ 00 ([7]), the function (/i/A)°'c<;i?(A, /i) remains a single-valued function of A^ and 
/i^, but it develops an essential singularity at A^ = 00, /z^ = 00. 
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In the present paper we shall consider sG model which requires putting together 
the two chiralities. To this end we shall need to consider not only the asymptotical 
region — )■ oo, but also — 0. Analysing the function ijjfi{X,fi) one concludes 
that the following limits exist 

oo 

(18) I \ira h* {{Car \) p. /3Leen(A) = E A"+^^-^/3:_,,. , 

oo 

\ \im cl{{Car\) ^ 7:crccn(A) = E A-"+2^7™,, • 

For the moment we do not know how to use these operators, but one thing is clear: 
they create highly non-local fields. 

The L-operator depends on ( and in a symmetric way. That is why another 
scaling limit is possible: 

(19) 2^*(A) = lim h^iCayX), 2j*{X) = lim c*((Ca)-'^A) , 

a— i>0 a->-0 

which allow the power series at A — t- 0: 

oo oo 

(20) ^*(A) :^ A'^^i-i, 7*(A) ^ >^'^%-i ■ 

The resulting operators act as follows 

^*(A) : V,®V,^2-i.(s-i) ^ V,®V„^2-i.s, 

These operators are obtainable from the previous ones by the substitution 

{a, A} ^ {2-a, A"^} . 



7*(A) : ^a®'^ a^2^{s+\) "aWi'a+aJ 



The proof goes through considering the scaling hmit (IT^ and studying the func- 
tion 

u}R{\,fi) = -( lim Wrat,n((Ca)"''A, (Ca)"» - Vuj{\/ fi, a)) . 

4 scaling 

The analysis is parallel to the one performed for the first chirality in IV. The func- 
tion {fi/ X)"u!ji{X, fi) is a single-valued function of A^, /i^ with essential singularities 
at A2 = 0, = 0^ 

The operators /3s(.j.een(A) and 7screen(A) are introduced similarly to the first chiral- 
ity. 

_ oo _ 

(21) I hm b*((Ca)-A) p. ^Leen(A) = E A"-^^^™,, , 

\- ^ CO '-' 
oo 

i hm c*((Ca)-A) p. 7s™(A) = E A^""" ^^7™,, • • 

^ a— >-0 A- -> oo 'J 

For both chiralities we have 

u}R{X,fx) wo(A//i,a), ujr{\,ij,) ujo{\/n,a) . 



So, the naive n — )• oo limit is reproduced. 
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3. Inhomogeneous six vertex model and sine-Gordon model 

We want to put two chiral models together and to make them interacting. Ac- 
cording to the previous discussion the lattice analogue of chiral CFT is the XXZ 
model. So, the two non-interacting chiral models correspond to the lattice contain- 
ing two non-interacting six vertex sublattices. As an example we shall consider the 
"even" or "odd" sublattices, consisting of lattice points with coordinates (j, m) such 
that j, m are both even or both odd. It is well-known [35] how to force these two to 
interact: one has to consider a six vertex model on the entire lattice with alternating 
inhomogeneneity parameters. We denote S = SUS, M = MUM and introduce 



rv 

oo 



2n 



S,M 



m=l 



Here 



^i,m(C) = Q 



Note that for = oo the inhomogeneous lattice reduces to two non-interacting 
homogeneous lattices. 

The functional Z^^^ is defined by 



^full|^2a5(0)Q 



TrsTrM(Ts,Mg'"^(°)0 



(22) 



We set K = here. The methods of II, III allow one to compute this functional with 
inhomogeneities in both directions. In particular, from II one concludes that the 
annihilation operators split into two parts: 

b(C) = b+(C) + b-(C), c(C) = c+(C) + c-(C) , 



b^(C) 

Let us take n = 
obtains: 

(23) 

We have 
where 



p=0 p=0 

oo. Then using the remark from the previous section and II, III one 



;7full r„2a5(0)(n\ ^ Trg (c"'"" {q 



2a5{0)(r)j^ 



n 



full 



^0 







(27ri) 




c^o(C/e,«)c^'(Ob1C)^^, 



and similarly for Vf2. The contour T± goes anticlockwise around Co"^- 
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As it has been said in the introduction, ideally we would like to start from a non- 
critical (XYZ or SOS) lattice model and to obtain the relativistic massive model by 
the usual scaling limit near the critical point. Since we do not have the necessary 
formulae to do that, we have recourse to the scaling limit of an inhomogeneous 
model by the procedure of [35]. Here we have an ideal situation from the point of 
view of QFT. Namely, we have the ultraviolet cutoff (lattice), the infrared cutoff 
(a finite number n of sites in the Matsubara direction), and the physical quantities 
(the values of the the functional Z^"^^ on quasi-local fields) are exactly computed 
with finite cutoffs. 

Like in the homogeneous case, let us introduce the step of the lattice a, and 
consider the scaling limit fll4p : n — )■ oo, a — )■ 0, na = 2itR fixed. We require further 
that Cq""^ — )■ 0, so that 

(24) M = 4a-\-^^'' fixed. 

The parameter M is a mass scale which has the meaning of the sG soliton mass [35] . 
The famous formula relating the soliton mass to the dimensional coupling constant 
H [36] reads in our notation as 



(25) ^ 



AC 



In this paper we consider a further limit i? — oo. In that case the sG partition 
function is obtained from Z^^^. 

Let us give some explanation at this point. The subject of study in [35] is the par- 
tition function of the sG model on the cylinder. There are two possible approaches 
to this partition function which correspond to two Hamiltonian pictures. In the first 
one, the space direction is considered as space and the Matsubara direction as time 
(space channel). One has scattering of particles and describes the partition function 
by the Thermodynamic Bethe Ansatz (TBA) [38]. This approach has an advantage 
of dealing with known particle spectrum and S'-matrices. It also has a disadvantage, 
because as usual in the thermodynamics one has to deal with the density matrix, 
which is a complicated object even in integrable cases. 

The paper [25] uses an alternative picture: the Matsubara direction is space, and 
the space direction is time. In this approach, the partition function is described by 
the maximal eigenvalue of the Hamiltonian of the periodic problem for the Matsub- 
ara direction (Matsubara channel). The advantage of this approach is clear: one 
deals with the pure ground state instead of the density matrix. The disadvantage is 
that describing eigenvalues in the finite volume is a difficult problem. This problem 
is addressed in [53] . 

More precisely, it is proposed in [53] to obtain the sG partition function as the 
scaling hmit of TrsTrM(7s,M)- As in the present paper, it is important to be able to 
control the computations starting from the lattice and from finite n. In our opinion, 
the main achievement of [35j is not in rewriting the Bethe Ansatz equation for 
the Matsubara transfer matrix in the form of a non-linear integral equation, but in 
extracting a main linear part and inverting it. The resulting Destri-DeVega equation 
(DDV) has several nice features. First, it allows the scaling limit and the mass scale 
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appears. Second, after the scaling limit the DDV equation clearly allows the 
large R expansion. Third, the scattering phase of the sG solitons appears in the 
DDV equation in the Matsubara channel. The last property allows the identification 
with the space channel. In particular, M happens to be equal to the mass of soliton. 

Now we come to the most important point of this paper. We wish to define the 
creation operators appropriate for taking the scaling limit to the sG theory. The 
issue is similar to the one in conformal perturbation theory, where one needs to 
prescribe a way how to extend the descendants in CFT to the perturbed case. It is 
claimed in [HI |23] that after subtracting the divergencies these operators are defined 
uniquely: possible finite counterterms can be dismissed for dimensional reasons, at 
least in the absence of resonances [THl [23] • The latter condition is satisfied in our 
case if u and a are generic. 

Let us define the operators hQ{() and Cq{() by the same formula as in the homo- 
geneous case (|T0|) . Starting from these operators we define the creation operators 

oo 

p=l 

and likewise for Cq(^). 

The operators b^* c^* create quasi-local fields. Notice that Z^'^ is defined on the 

quotient space Wquo because the t*-descendants do not contribute to it. 

There are two sorts of chiral operators, one living on the even sublattice and the 
other on the odd sublattice. We have to make some combinations which will give 
finite answers for the interacting model. At the same time, we want this combination 
to correspond to our intuitive idea that we have to subtract perturbative series. Let 
us explain that the correct combinations are given by the following. 

(26) b+*(C) = e-^^br(C)e^o' , c+*(C) = e-^^cr(C)e™'" , 

b-*(C) = e-^o"bo*(C)e^«" , c-*(C) = e-^o-Co*(C)e^°" • 

Corresponding creation operators which create the quasi-local operators are defined 
by 

oo 

p=l 

In what follows we shall be interested only in the case of an equal number of b"^* 
and c+*. Let us examine how the descendants of this form depend on (o, taking the 
simplest case 

br 4* K* cr (g^-^W) = e-^^"-^o--b+* c+; h,; Co;:(g^"^(°)) . 

It is easy to see from the definition in II that the second factor in the right hand 
side is a rational function of Co regular at oo, 

(27) b+* c+; bo-: Co;:(g^""(°)) = (Oo + Co-'Oi + Co-'02 + ■ ■ •) q"''^'^ • 

The rationality holds true after application of e"^o ^~^o . The expansion (^7^ looks 
as a perturbative series for the action ([1]), because in the scaling limit (q'^ comes 
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accompanied by a and Co ^'^ has the dimension of (mass)^'' i.e. that of fi^ 
(see flM|) ). Notice that the property fl27j) would be spoiled if we apply e~^o -^o^ 
as well, because it will pick up ujo{(, a) near ('^ = (^"^^ and the function uo{(, a) has 
asymptotics at ^ — )■ cx3 containing both (^""2™ and C ~- 

On the descendants created by fl26|) . the value of Z^^" remains finite in the scaling 
limit. So we conclude that they create renormalised local operators. We cannot 
say anything about the finite renormalisation, but it can be taken care of by a 
dimensional consideration (see below). 

We conjecture that the following limits exist 

(28) |limb+*(C) „^ /3*(K) + ^™(C/m), 

scaling (, — > oo 

i hm C+*(C) r(K)+7:creen(C/M), 
scaling f — !> oo 

i hm b-*(C) ^ r(C/A^)+/3:creen(MC), 

scaling (, U 

1 hm C-*(C) ^ r(C/M)+7:creen(K), 
^ scaling C 

where by limscaiing the scaling f|T^ . fl24l) and i? — > oo are implied. 

In these formulae we denote by the same letters the operators in the sG model as 
they were denoted in the CFT. We shall not consider the screening operators,while 
we use the coefficients /32j-i so forth, defined as in f|T71) and fl20|) . to consider 
the descendants 

(29) ^;-+7/-/3J+7/-$.(0). 

This is the field in the interacting model which goes to corresponding descendant in 
the conformal limit — ^ , and which does not develop finite counterterms. The 
latter are forbidden by dimensional consideration. So, this is exactly the definition 
which we were supposed to use from the very beginning. Notice that the appearance 
of fj, in the formulae (|28|) is a consequence of consistency with the conformal limit 
due to (1251). 

Now it is easy to compute the normalised vacuum expectation value of the de- 
scendant fl29l) for the sG model. It is obtained by taking the scaling limit of Z^^^ 
and computing the asymptotics of Uo{(,a) for ^ oo and C which is done 
simply by summing up the residues in appropriate half-planes: 

uoiC, a) ~ - V cot — (2n - 1 + z/a) + i V C""^" tan — (a - 2n) . 

V ^-^ 2u ^-^ 2 

n>l n>l 

The first part corresponds to the expectation values of operators created by the 
coefficients of the expansions (fT7|) . (!20|) . These are the expectation values in which 
we are interested in this paper. The second part of the asymptotics corresponds to 
the operators (fT8|) . (!2T|) At present we do not know the meaning of their expectation 
values, but we hope to return to them in the future. 
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Introducing the multi-indices as described in the introduction we obtain: 

($a(0))sG \l^J 

X ^Id^+I+I^-I) JJ cotf(^ + «) n cotf(^-«). 

2n-le/+ 2n-le/- 

Now we have to recall the definition in IV: 

(31) /3L-l = ^2n-l(a)/32'ri% 72n-l = ^2n-l(2-a)7S^':T^ 

where 

B (r.\ ^ r^-i r(f + ^(2n-i)) 



IV in 



-l)!r(f +(^(2n-l)) 



with 



Similarly we get for the second chirality 

The main formula ([5]) follows immediately. 

Before concluding this paper, let us say that in principle our approach can be 
applied to the computation of one-point functions of descendants for finite radius in 
the Matsubara direction (finite temperature, in other words). However, this would 
require a detailed study of the DDV equation and the equations for the function ur 
for the sG model in finite volume. 
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